Vortex sheet dynamics and turbulence 
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The nonlinear evolution of a vortex sheet driven by the Kelvin-Helmholtz instability is characterized by the 
formation of a spiral possessing complex stretching and intensity patterns. We show that the power energy 
spectrum of a single two-dimensional vortex sheet tends to the usual fluid turbulent spectrum, with an exponent 
of -3. Using numerical simulations and asymptotic methods, we demonstrate the relation between this power 
law and the singularities in the geometry and vorticity distribution of the sheet. 
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In 1941, Kolmogorov formulated a theory of fluid turbu- 
lence, based on the idea of statistical similarity (H. In par- 
allel, Onsager proposed that a singular enough velocity field 
could explain dissipation without the "final assistance of vis- 
cosity" yl (A thorough review on Kolmogorov theory can 
be found in the book by U. FrischM. Onsager's work on tur- 
bulence is accounted for in Ref. |5j.) These two complemen- 
tary approaches, the one statistical, the other dynamical, can 
be taken as constitutive of further theoretical developments on 
turbulence. One difficult problem, whose solution would con- 
tribute to understanding the physics of turbulent flows, is to 
establish the relation existing between dynamics and statistics, 
in particular flows. Following the theoretical ideas first intro- 
duced by Townsend |6], Lundgren proposed a relation of this 
type in his model of intermittent structures at high Reynolds 
numbers |H ■ Lundgren demonstrated that solutions of the 
Navier-Stokes equation in the form of slender spiral vortices 
yield Kolmogorov spectrum, through a mechanism of enstro- 
phy intensification and creation of fine scales by differential 
rotation and stretching. In two dimensional flows, Gilbert [9] 
studied the wind-up of a vortex patch by a strong point vortex, 
and found that an energy spectrum E(k) ~ k~ y with 3 < y < 4, 
is obtained from the fractal structure of the spiral separating 
potential and vortical regions. The goal of this Letter is to 
study the relation between geometry, vorticity distribution and 
spectral properties of the flow generated by the winding-up of 
a shear layer subject to the Kelvin-Helmholtz instability. 

In two dimensions, and in the limit of vanishing viscosity, a 
shear layer becomes a line of velocity discontinuity. The Euler 
flow associated with a vortex sheet reduces to the Hamiltonian 
dynamics given by the Birkhoff-Rott equation 

jZ(s, t) = ^- jds'K(s\ t)K(z(s, t) - Z(s', f)) . (1) 

In Eq. the vortex sheet is represented by the parametric ex- 
pression of the curve z = Z(s, t), where z — x + iy is a complex 
coordinate and s is the arclength, ds = \dZ\. The Biot-Savart 
kernel is defined by K(z) — l/z or K(z) = cot(7rz)/2 in the 
case of periodic boundary conditions in the x-direction. In 
J!) k — F s = dF/ds (subscript s denotes derivative) is the 
intensity of the sheet, where the circulation F = T(s, f) is a La- 
grangian variable satisfying dF/dt — 0; § denotes the Cauchy 



principal value integral, and the over bar denotes the complex 
conjugate. A sinusoidal perturbation of a straight vortex sheet 
evolves towards a singularity in finite time, as demonstrated 
by Moore in 1979 fvk If3| . probably related to the incipient 
rolling-up of the sheet into a spiral. In order to integrate Q 
in time, it is then necessary to regularize the kernel. A simple 
method consists in cutting the vortex interaction at a distance 
5, therefore replacing K ~ l/z by K s = z/(\z\ 2 + S 2 ) fil. 
Solutions of the regularized Birkhoff-Rott equation, with K 
replaced by Kg, tend to solutions of the Euler equation in the 
limit 6 — > 1 15]. It is worth noting that the dynamics of the 
desingularized Q remains Hamiltonian. As a consequence, 
this model is particularly well suited to investigate the dy- 
namics of the turbulent cascade in one of the most basic flow 
structures, the Kelvin-Helmholtz unstable vortex sheet. 

We used the blob vortex method of Krasny fl4ll to compute 
the dynamics of a shear layer in a periodic domain of length 
L. We take units such that L = 1 and the circulation is in 
the interval F e [0, 1] (or the unperturbed velocity disconti- 
nuity equal to one). The initial perturbation, as for example 
the superposition of small amplitude modes, evolves towards 
the formation of a system of spirals winding around the main 
one. These structures are driven by a cascade of secondary 
Kelvin-Helmholtz instabilities, as shown in Ref. |16]. Fig- 
urenp resen ts the vortex sheet formed from an initial condi- 
tion of two modes with wavelengths 1 and 1/4, together with 
the strain rate S{s, t) and the intensity k(s, t). The strain rate is 
defined by S = (l/As)dAs/dt, where As — > 0, is the arclength 
separation of two neighboring points on the vortex sheet. Di- 
minishing the value of 5, the relevant physical limit is 5 — » 0, 
results in a drastic increasing of small scale features of the 
sheet winding 1 17] (see also Il8l0 . We study numerically the 
generation of the small scales resulting from the short wave in- 
stabilities, in the presence of nontrivial stretching and vorticity 
concentration effects. Indeed, arbitrarily small scale perturba- 
tions might grow, depending for instance on the local strain 
rate 11611 . and regions of strong localized vortex intensity de- 
velop as observed in Fig.^b). The origin of these small scale 
features, can be traced back to the highly oscillatory behavior 
of the vortex sheet stretching, whose dynamics is essentially 
Lagrangian and nonlocal (Fig.Qc)). The presence of strong 
velocity gradients (in the limit 5 — > 0) possessing a wide range 
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FIG. 1: (Color online.) Spiral formed at the nonlinear stage of the Kelvin-Helmholtz instability (N = 16384, 5 = 0.024, t = 0.93). (a) 
Snapshot of the vortex sheet; (b) the intensity k(s, t) concentrates in localized regions; (c) the strain rate S(s,t) (in blue) shows a complex 
structure of alternating compression and stretching zones (in red, the vortex sheet). 



of space and time scales, responsible for the appearance of an 
increasingly complex velocity field, can be considered to be 
the root of a turbulent flow. We found indeed, that the flow 
created by this single shear layer becomes turbulent, in the 
sense that it is characterized by a well defined power law en- 
ergy spectrum Ek ~ k~ y (for k <K 2n/6, and 3 < y < 4), in 
agreement with full Navier-Stokes simulations with random 
initial vorticity field fFolEoll . 

It is then natural to look at the energy spectrum, E% = E(k) 
of the sheet. However, in order to avoid the arbitrariness in the 
origin of coordinates we define a one-dimensional spectrum 
from the y-direction average of the velocity v(x,y) correlation 
function c(g) , 

l 

E k = E(k) = intctf)] = \ J <%c(?)e- ik t , (2) 

o 

where Tk is the Fourier transform, and 
i l 

= Yi J d y J dxv ( x + &)0 • v ( x > y) ■ ( 3 ) 

-/ 

(The averaging length I is chosen munch larger than the spiral 
size.) This definition, similar to the one proposed by Saffman 
and others 0, 0, E2I1 . satisfies the condition that the con- 
served total energy E, is given by J Q dkE(k) = J. dk\Vk\ 2 = 
E, where, - Tk\v\. We measured the spectrum of the 
vortex sheet shown in Fig. ^ an d found that for k «: 2n/6, 
it is well fitted by a power law Ek ~ AT 2 94 close to the law 
observed in developed decaying two-dimensional turbulence 
In fact, the actual value of the exponent follows the evo- 
lution of the initially regular vortex sheet. The value y » 3 is 
reached after the development of the secondary structures (see 
the insert in Fig.^. At the early stages of the winding, inter- 
mediate values 3 < y < 4 are measured, with a low wavenum- 
ber range dominated by y « 3, and a large wavenumber range 
having a steeper spectrum. The y w 3 range of k increases 
with time, and at t — 0.386 it extends almost over all the 
available range (0 < k < 2n/6). In order to obtain a sig- 
nificant range in the turbulent spectrum we proceeded for the 
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FIG. 2: (Color online.) The energy spectrum J2j of the flow associ- 
ated with the vortex sheet shown in the insert (N = 2 15 , S = 0.01, 
/ = 0.386, quadruple precision computation). 



first time to a quadruple precision computation (33 significant 
digits), using 6 = 0.01 and = 2 15 vortices. 

As it is well known, the superposition of velocity disconti- 
nuities yields a spectrum in Ek ~ k~ 2 |6j, while a superposi- 
tion of vorticity discontinuities yields Ek ~ AT 4 12 ill , the clas- 
sical enstrophy cascade picture gives Ek ~ k~ 3 12311 . It would 
be interesting, as suggested by Moffatt 1 24], to establish a rela- 
tionship between the fractal structure of the rolled shear layer 
and the observed energy spectrum, characterized by an inter- 
mediate value of the exponent. (It is worth mentioning that for 
finite 5 the velocity field is continuous.) The fractal dimension 
df of the spiral in Fig.^is found to be df- 1.44. This value 
can be compared to d? = 1.33 measured in Ref. 1 25], for a 
single Kelvin-Helmholtz spiral and for a large value of the 
smoothing parameter 5 = 0.25. Assuming that the formula 
y = -5 + d F relating the energy exponent to the fractal di- 
mension, derived by Gilbert 0] in his study of a spiral vortex 
patch, holds in the present case, we obtain a value y as 3.56, 
somewhat larger than the observed one. The difference may 
be attributed to the nonuniformity in the distribution of the 
vorticity. 

It is then tempting to conjecture that the characteristic ex- 
ponent of the spectrum y reflects, not only the fractal spiral 
structure of the sheets, but also the existence of singularities 
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FIG. 3: (Color online.) Shape and spectrum of the vortex sheet, (a) The spiral shown in Fig.Qafter periodization in the x direction; (b) power 
spectrum of Z(s), with exponent 2a + 2 = 2.45; (c) power spectrum of k(s), with exponent 2{S + 2 = 0.66. 



in both the shape and the intensity of the vortex sheet. Remark 
that the existence of such type of singularities is not related to 
eventual singularities in Euler flows, the velocity field outside 
the sheet is always regular in our case. Indeed, we can use the 
Biot-Savart equation for the velocity field to relate y with the 
characteristic exponents a and p of the shape Z(s) ~ s a and 
intensity k(s) ~ singularities, respectively: 



I"'! 



dx- 



£-GS+l)/a _ 



(4) 



Therefore, Eu = | v^| 2 ~ k~ y implies y = 2(J3 + l)/a. The 
exponents a and p can be directly measured from the data of 
Fig. [5] where we display the shape Z(s) (after subtraction of 
T(s), in order to obtain a periodic function Z(F + 1) = Z(T), in 
terms of the parameter F) and its spectrum together with the 
one of the intensity. In Fourier space we have |Zi| 2 ~ k~ 2a ~ 2 
and \K k \ 2 ~ k- 1 ?- 1 , it follows that a = 0.22 and p = -0.66, 
giving y — 3 in agreement with the exponent measured from 
the energy. 

We found that the intensity concentrations shown in Fig.^ 
are related with regions characterized by the strong advec- 
tion of secondary structures. Indeed, during the evolution of 
secondary instabilities, filamentary sheets are advected by the 
previous formed intense vortices (the central tightly wounded 
region of the spirals). A simple model describing this situa- 
tion, consists of a fixed strong point vortex advecting an ini- 
tially straight weak intensity vortex sheet. 

The evolution of the sheet z = Z(s, t) of intensity k(s, t) - 
dT(s,f)/ds, s is the arclength, is given by a Birkhoff-Rott 
equation with an external velocity field due to the point vortex 
of circulation Fq. 



— Z(s,t) = 

dt 2niZ(s,t) 



1 £ K(s',t)ds' 
2m J Z(s,t) -Z(s'j) ' 



Initially one may consider that the sheet is located at a distance 
d in the y direction from the central point vortex, Z(s, 0) = 
s + id. If one neglects the self-interaction, the advection of the 
sheet is described by the curve, Z = za(s, t) — (s + k/)e' n( ' 5) ', 
where Q(s) = Tal2n{d 2 + s 2 ). However, a singularity can- 
not result from the sole advection, and the effect of self- 
interaction is needed. Direct integration of actually shows 



that, in the limit of 6 — > 0, the curvature diverges at a point 
So, where it changes sign, leading to the eventual formation 
of a new spiral (Fig.0}. In order to identify the nature of this 
singularity we assume that the behavior of the vortex sheet 
around s = sq, can be analyzed within a local approximation 
Il26ll . therefore neglecting the contour integral in the computa- 
tion of the principal value of the Birkhoff-Rott term, 



d - i a 

—Z(S, f) = —=7—T + , 

at Z(s, t) Z s (s, t) 



(6) 



where we introduced non-dimensional variables based on the 
time unit 2nd 2 /Fo and the space unit d. The non-dimensional 
parameter a = nKod/To <k 1 measures the circulation ratio be- 
tween the central point vortex and the sheet. In the following 
we put the point z = Z(sq, to) corresponding to the singularity, 
at the origin sq = and fo = 0. Moreover, the transformation 
Z(s, t) = r(s, t) exp[i(/i(s, t) + if2(s)f], allows to separate the dy- 
namics of the complex amplitude re 1 * from the pure rotation 
at frequency £2, and leads, after substitution into (jSJl, to the 
two real equations, 



— r - 



d 



■a,t) 



2 ' 



-O + -r \ + a 



r 2 + r 2 ((p s + Q. s f) 2 



(7) 



(8) 



We assume now, that the form of the sheet near de singu- 
larity can be written as r = ro + sgn(st)R)(s, t) and cp = 
<f>o + H(t)<S>)(s, t), with the regular parts determined by the ad- 
vection r Q (s) = |z | = [1 + (s + s ) 2 ] 1/2 and <f> (s) = arg(z ) = 
arctan(s + sq) + n/2 + £2fo, and where the singular parts R(s, t) 
and 0(s, t), should possess strong derivatives near s * 0. This 
can be satisfied if R ~ \s\ q , with the condition q < 1. We in- 
troduced the sign sgn(sf) and the Heaviside H(t) functions, to 
reproduce the jump in the curvature |AC(f)l, as observed in the 
numerical computations of Fig. |4] together with the smooth- 
ness of the phase <D>(i, f) before the singularity. Substituting 
these expressions into J7I8L and neglecting small terms in 
\s\ <k 1] (we can a posteriori verify the consistency of these 
assumptions), we get, 

8 sgn(jp d sgn(j) O s + Cl s t 
—R-a , — <D = — R + a , (9) 

dt R s dt r\ R 2 
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of strong secondary singularities, associated with jumps in 
the curvature, by the simultaneous action of the advection and 
self-interaction of the vortex sheet. A simple model describes 
the structure of the sheet near the singularity in the form of a 
square root cusp Z(s) ~ sgn(s)|s| ? , with an exponent q — 1/2. 
This value is sufficiently small to be in agreement with the 
requirement that 3 < y < 4. If one assumes that q ~ a, the 
condition (J3+ 1)/2 < q ~ a < 2/3 should hold. The formation 
of singular concentrations of vorticity in regions with diverg- 
ing curvature, constitute a dynamical mechanism towards a 
turbulent dissipation in the limit of vanishing viscosity. 



FIG. 4: (Color online.) The maximal curvature jump as a function of 
time for various values of 5 e (0.013, 0.005), showing the approach 
to a singularity in the limit 6 — » 0, during the advection of the sheet, 
from numerical solutions of Inserts: (top) Snapshot of the vortex 
sheet with 6 = 0.005 at / = 1.8; and (middle) the corresponding 
vortex intensity k as a function of the arclength s. 



whose solutions, singular at s ~ 0, are, 
R = (4a|sf|) I/2 , 



and 



n,|j|log(l/|a1)-^^(4flM) 1/2 



(10) 



t + f(st). (11) 



where f(st) is an arbitrary smooth function, that can be ne- 
glected. We note that the solution of l|9} naturally leads to 
q = 1/2, and to the scalings: R s » ro<f> s for t <k 1, justifying 
the assumption that the second term in the denominator of (fTJ- 
|HJ is negligible near the singularity. It is also worth noting that 
this singularity is stronger than the one found by Moore (3, 
whose characteristic exponent is q — 3/2 > 1, the difference 
coming from the presence of the advection velocity field. 

In this Letter the k~ 3 energy spectrum of two dimensional 
turbulence is recovered in a simple model of a single, regu- 
larized, vortex sheet evolving according to a Hamiltonian dy- 
namics. The vortex sheet is subject to a cascade of Kelvin- 
Helmholtz instabilities, leading to spiral windings at different 
scales. The flow is characterized by an intermittent distribu- 
tion of the vorticity intensity along the vortex sheet, and by a 
complex sequence of compression and stretching regions re- 
sulting from the nonlocal Lagrangian dynamics of the strain. 
The sole consideration of the fractal structure of the spiral 
windings is not enough to account for the observed spectrum. 
A simple dimensional argument using the Biot-Savart equa- 
tion, allows to relate the exponents of the singularities in the 
shape and the intensity of the vortex sheet, with the spec- 
trum exponent. We propose a mechanism for the formation 
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